Draft version January 13, 2013 

Preprint typeset using MTgX style emulateapj v. 05/04/06 



A CHANDRA STUDY OF THE LARGE-SCALE SHOCK AND COOL FILAMENTS IN HYDRA A: EVIDENCE FOR 
SUBSTANTIAL GAS DREDGE-UP BY THE CENTRAL OUTBURST 



Myriam Gitti 



and Michael W. Wise 



o 

(N 

X> 

<D 

(N 

o 

u 

43 

O-i 

6 



> 

m 
in 

o 



PaulE. J. Nulsen , Laurence P. David , Brian R. McNamara 

Draft version January 13, 2013 

ABSTRACT 

We present the results of a Chandra study of the Hydra A galaxy cluster, where a powerful AGN outburst 
created a large-scale cocoon shock. We investigated possible azimuthal variations in shock strength and shape, 
finding indications for a weak shock with a Mach number in the range ~1.2-1.3. We measured the temperature 
change across the shock front. However, the detection of a temperature rise in the regions immediately inside 
of the front is complicated by the underlying temperature profile of the cluster atmosphere. We measured the 
global temperature profile of the cluster up to 700 kpc, which represents the farthest measurement obtained 
with Chandra for this cluster. A "plateau" in the temperature profile in the range ^70-150 kpc indicates the 
presence of cool gas, which is likely the result of uplift of material by the AGN outburst. After masking the 
cool filaments visible in the hardness ratio map, the plateau disappears and the temperature profile recovers 
a typical shape with a peak around 190 kpc, just inside the shock front. However, it is unlikely that such 
a temperature feature is produced by the shock as it is consistent with the general shape of the temperature 
profiles observed for relaxed galaxy clusters. We studied the spectral properties of the cool filaments finding 
evidence that ~1O U M0 of low-entropy material has been dredged up by the rising lobes from the central 30 
kpc to the observed current position of 75-150 kpc. The energy required to lift the cool gas is >2.2 x 10 60 
erg, which is comparable to the work required to inflate the cavities and is ^25% of the total energy of the 
large-scale shock. Our results show that the AGN feedback in Hydra A is acting not only by directly heating 
the gas, but also by removing a substantial amount of potential fuel for the SMBH. 

Subject headings: galaxies:clusters:general - galaxies: clusters: individual (Hydra A) - cooling flows - inter- 
galactic medium - galaxies:active - X-rays: galaxies: clusters 



1. INTRODUCTION 

Observational and theoretical evidence has been growing 
in the past decade in favor of the existence in galaxy clus- 
ters of a feedback mechanism that prevents cool cores from 
establishing "cooling flows" at the rates predicted by earlier 
X-ray observations (e.g., Peterson & Fabian 2006 for a re- 
view). The dominant cD galaxies, which are present at the 
cluster center in all cool core clusters, host the most massive 
black holes in the local Universe and usually show nuclear 
activity. They accordingly provide a natural feedback mech- 
anism for the regulation of the cooling process. Feedback is 
also required to suppress the overproduction of massive galax- 
ies predicted by dark-matter-only simulations and to break the 
self-similarity of clusters (e.g., Benson et al. 2003). The na- 
ture of this feedback, vital to our understanding of galaxy and 
structure evolution, is one of the most important unresolved 
questions in extragalactic astronomy. Based on the detec- 
tion of cavities and AGN-driven shocks, the primary source 
of feedback in cluster has been identified as radio galaxies 
acting through outbursts and accompanying energy injection, 
likely intermittent, from the central AGN (e.g., for a review, 
McNamara & Nulsen 2007 and references therein). In galaxy 
clusters, where cooling rates should be highest, the current 
generation of X-ray observatories Chandra and XMM-Newton 
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have shown that there is not a significant amount of gas cool- 
ing below about one third of its virial temperature (Peterson et 
al. 2003; Kaastra et al. 2004). Images from these telescopes 
also reveal highly disturbed structures in the cores of many 
clusters, including shocks, cavities and sharp density disconti- 
nuities. At radio wavelengths, it is clear that AGN jets are the 
cause of many of these disturbances. The incidence and vari- 
ety of bubbles, cavities, shocks, and ripples observed both in 
the radio and in X-rays in galaxy clusters provides direct evi- 
dence of the widespread presence of AGN-driven phenomena 
(e.g., Fabian et al. 2003 for a discussion of the Perseus clus- 
ter properties in terms of AGN-generated viscously-damped 
sound waves; Nulsen et al. 2005 for shock heating in Hy- 
dra A; Gitti et al. 2007 for a study of giant cavities created 
by the most powerful AGN outburst currently known; Birzan 
et al. 2008 and Diehl et al. 2008 for a survey of cavities 
and the implied cavity heating rates). Such AGN feedback 
has a wide range of impacts, from the formation of galax- 
ies, through to the explanation of the observed M^-cr relation 
(which indicates a causal connection or feedback mechanism 
between the formation of bulges and their central black holes, 
e.g., Magorrian et al. 1998), to the regulation of cool cores. In 
most cases, the energy introduced by the AGN is more than 
sufficient to counteract putative cooling flows (Birzan et al. 
2004, 2008; Rafferty et al. 2006). However, the details of 
how the feedback loop operates are still unknown. Only by 
studying striking examples of interaction between the central 
radio galaxy and the ICM can we understand why cooling and 
star formation still proceeds at a reduced rate, and this is likely 
to reveal the coupling between AGN feedback and the ICM. 

The galaxy cluster Hydra A has a well-known, large-scale 
system of X-ray cavities embedded in a "cocoon" shock sur- 
rounding the central, powerful radio source (McNamara et al. 
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TABLE 1 
Shock properties 



Sector Shock radius Radius variation Density jump Mach number 



01-02 (°) 


Re (") 


R B] - R e 2 (") 




M 


0-30 


208 


201-222 


1.33 


1.22 


30-60 


257 


222-303 


1.45 


1.30 


60-90 


332 


314-349 


1.33 


1.23 


90-120 


371 


364-342 


1.38 


1.26 


120-150 


301 


342-271 


1.46 


1.32 


150-180 


260 


264-265 


1.40 


1.27 


180-270 


285 


260-253 


1.15 


1.10 


270-360 


205 


233-202 


1.35 


1.23 



NOTE. — Results of the fits of a broken power-law density model to the surface bright- 
ness profiles extracted along different sectors. The sector aperture (from RA. 9\ to #?) is 
indicated in the first column. The shock radius at the mid-angle of the sector (R§) is indi- 
cated in the second column, whereas the shock radii at the starting (R$ i ) and ending (Ro [ ) 
angles are indicated in the third column. Such radial distances are measured from the cluster 
center. The corresponding density jump and Mach number (Ai) are shown in the forth and 
fifth columns, respectively. Best-fit statistical errors are on average ~5%. 

2000, Nulsen et al. 2005). It is considered one of the proto- 
types of cool core clusters with cavities, which has served as 
an early test of the AGN feedback paradigm, and it has been 
extensively studied both in the radio and X-rays (Taylor et al. 
1990; McNamara et al. 2000; David et al. 2001; Nulsen et al. 
2002, 2005; Lane et al. 2004; Wise et al. 2007; Simionescu 
et al. 2009a, 2009b, Kirkpatrick et al. 2009). By analyzing 
the archival —200 ks Chandra exposure in this paper we study 
the azimuthal properties of the large-scale shock and attempt 
to measure the temperature jump associated with the shock 
in different directions. We also perform a detailed spectral 
analysis of the cool X-ray filaments extending out to 150 kpc, 
finding evidence for extensive mass dredge-up from the cen- 
tral 30 kpc, which affects the global temperature profile of the 
cluster. 

With H = 70 km s" 1 Mpc" 1 , and il M = 1 - fL A = 0.3, the 
luminosity distance to Hydra A (z=0.0538) is 240 Mpc and 1 
arcsec corresponds to 1.05 kpc in the rest frame of the cluster. 

2. CHANDRA OBSERVATIONS AND DATA REDUCTION 

Hydra A has been imaged four times by Chandra ACIS for 
a total exposure of 240 ks. The two shorter exposures (Ob- 
sIDs 575 and 576, which collectively comprise only — 17% 
of the total exposure) were taken early in the Chandra mis- 
sion when the ACIS detector was operated at the higher focal 
plane temperature of -1 10°C. Due to the higher quality cali- 
bration at -120°C, the spectral analysis presented below was 
performed on the two longer, more recent exposures: —97 ks 
collected on 2004 January 13 (ObsID 4969) and -99 ks col- 
lected on 2004 October 22 (ObsID 4970) with ACIS-S. We 
use data from the S3 and S2 CCDs to study the central part 
of the cluster emission where the cavity system and the radio 
source are located, and from the 12 and 13 CCDs to measure 
the temperature and surface brightness in the cluster outskirts. 
Each dataset was individually reprocessed with CIAO ver- 
sion 4.1 using CALDB 4.1.0 and corrected for known time- 
dependent gain and charge transfer inefficiency problems fol- 
lowing techniques similar to those described in the Chandra 
analysis threads 6 . Screening of the event files was also applied 
to filter out strong background flares. Blank-sky background 
files, filtered in the same manner as in the Hydra A image 
and normalized to the count rate of the source image in the 

6 http://cxc.harvard.edu/ciao/threads/index.html 




Fig. 1. — Residual map of the beta model subtracted surface bright- 
ness image from Nulsen et al. (2005). The green contours outline the 
330 MHz radio emission from Lane et al. (2004). Overlaid in white 
are the sectors considered to investigate the azimuthal variations in 
shock properties by means of the surface brightness profile. In each 
sector, the red curve indicates the position of the shock front deter- 
mined by fitting a broken power-law density model to the surface 
brightness profile (see Fig. |2]and Sect. 13. 1 I for details). In the sector 
180-270, where the surface brightness discontinuity is least evident, 
we have indicated the shock front with a dashed curve. 

10-12 keV band, were used for background subtraction. We 
identified and removed the point sources using the CIAO task 
WAVDETECT, with the detection threshold set to the default 
value of 10~ 6 . The final, combined exposure time for the two 
dataset is 174.2 ks. 

3. THE LARGE-SCALE SHOCK 

The observed feature interpreted as a shock front in the 
X-ray surface brightness of Hydra A surrounds the low- 
frequency radio lobes, and the correspondence between their 
shapes supports the interpretation of a cocoon shock of the 
radio source. Nulsen et al. (2005) estimated an age of the 
outburst of t s = 1 .4 x 10 8 yr and a total energy of E s = 9 x 10 60 
erg. The shock front is clearly aspherical with a complicated 
3-D shape affected by projection effetcs. In particular, the 
Hydra A radio source and cavity system is inclined at about 
40° to the plane of the sky, with the northern side lying closer 
to us (Taylor 1996, Lane et al. 2004, Wise et al. 2007). The 
cavity and shock geometry as well as the outburst history and 
energetics have been studied in detail in Wise et al. (2007) 
and Simionescu et al. (2009b). 

3.1. Azimuthal Variations in Shock Strength 

We have investigated the azimuthal variations in shock 
strength and shape by studying surface brightness profiles in 
different directions. In particular, we have divided the north- 
em cluster semicircle into six sectors of 30 degrees each, and 
the southern cluster semicircle, where the shock is less evident 
in the image, into two quadrants of 90 degrees each. Start- 
ing from the West with position angle (PA.) = 0° and count- 
ing counterclockwise, the sectors are labeled as: 0-30, 30- 
60, 60-90, 90-120, 120-150, 150-180, 180-270, 270-360 (see 
Fig. [TJ. We extracted the background-subtracted, exposure- 
corrected surface brightness profile along each sector in the 
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FlG. 2. — Background-subtracted, expos ure-correc ted Chandra surface brightness profiles extracted along the sectors outlined in Fig.0 The surface brightness is in units of counts 
cm -2 s -1 , with errors at 1 a. The radial axis shows the distance from the center of curvature of the annuli used to extract the surface brightness profile. Radial error bars show the bin 
sizes. The smooth curves show fits of the broken power-law density model which give the density jumps and corresponding Mach numbers reported in Table fTl 
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energy range 0.5-2.0 keV. The center of the annular regions 
used to extract each profile was chosen in order to best match 
the curvature of the radial bins with the shape of the shock 
front as seen in the image (Fig. [TJ. We then performed fits of 
a broken power-law density model to each surface brightness 
profile. This analysis assumes that the radius of curvature in 
the plane of the sky is the same as that parallel to the line- 
of-sight. A hydrodynamic model for the shock was made by 
initiating an explosion at the center of a hydrostatic, isother- 
mal atmosphere with a power-law density profile (see Nulsen 
et al. 2005 for more details on the shock model). The power- 
law index for the density profile of the unshocked gas was de- 
termined from the broken power-law fit to the surface bright- 
ness profile. The surface brightness profiles and the best-fits 
of broken power-law density models in the various sectors are 
shown in Fig. [2] where the radial axes indicate the distance 
from the center of curvature of the annular regions used to 
extract the profiles. However, in the following discussion we 
indicate the position of the shock in each sector (referred to 
as the "shock radius") in terms of the distance of the front at 
mid-aperture of the sector from the cluster center. The best -fit 
model for each sector is summarized in Table Q] and the cor- 
responding shock front is shown by red arcs in Fig. Q] We 
also investigated the effects of varying the radial binning and 
center of curvature of the annular regions used to extract the 
profile. We found that the fit results (i.e., shock radius and 
Mach number, .A/I) do not depend strongly on the particular 
choice of the extraction region nor on the radial range of the 
fit, with systematic variations <5% in each sector. 

The shock is clearly visible as a surface brightness jump 
in all sectors but the 180-270 sector (i.e. the SE quadrant), 
with a radius varying between ^205" in the E-W direction 
to ^365" in the N-S direction. The Mach number varies be- 
tween 1.20 and 1.32 in the sectors where the density jump 
is detected, whereas to the SE (sector 180-270) where the 
surface brightness discontinuity is less evident we estimate 
M. = 1.10. These results are consistent with Nulsen et al. 
(2005), and partially in agreement with Simionescu et al. 
(2009b), who report a detection of the shock front also in the 
S direction. 

3.2. Temperature Variations Across the Shock Front 

Our shock model predicts the emission-weighted temper- 
ature to rise across the front by ~8%, 10%, 15%, 20% for 
Mach numbers of 1.18, 1.23, 1.33, 1.40, respectively, reach- 
ing its peak at a distance of ^5-10% of the shock radius be- 
hind the shock and declining below the unshocked tempera- 
ture values inside radial distances of ^25% of the shock ra- 
dius. We used this information to optimize our selection of the 
pre-shock and post-shock regions in each sector and then ex- 
tracted the spectra in these regions using the SPECEXTRACT 
task, which also computes the corresponding event-weighted 
response matrices. Spectral fitting to a single absorbed apec 
model was performed in XSPEC version 12.3.1 in the 0.5- 
8.0 keV energy range. Abundances were measured relative 
to the abundances of Anders & Grevesse (1989) and a galac- 
tic hydrogen column of 4.68 x 10 20 cm" 2 (Dickey & Lockman 
1990) was assumed. As a general result, the post-shock re- 
gions are found to be hotter than the corresponding pre-shock 
regions, although due to the large error bars the temperature 
of the pre-shock and post-shock gas is consistent. Some ex- 
amples of these temperature measurements are presented in 
Sect. 6. 1 . A further attempt to bin sectors together in order to 
increase the statistics did not lead to any clear detections of a 




Fig. 3. — Chandra 0.5-7.0 keV mosaic of the Hydra A field from 
Wise et al. (2007). Overlaid in green are the annuli used to measure 
the azimuthally-averaged temperature profile shown in the left panel 
of Fig. [4] The blue sectors indicate the regions where the cool fil- 
ament is present. These sectors have been excluded in determining 
the temperature profile shown in the right panel of Fig. |4](see Sect. [5] 
for details.). The red contours outlining the 330 MHz radio emission 
from Lane et al. (2004) are also shown for comparison. 

temperature jump across the shock front (see Sect. 6. 1). How- 
ever, previous measurements (David et al. 2001; Simionescu 
et al. 2009b) indicate that the underlying global temperature 
profile has a peak close to the inner edge of the shock front, 
thus complicating the detection of a temperature rise due to 
the shock itself. To investigate this in more detail, we study 
the properties of the azimuthally-averaged cluster temperature 
profile in Sect. [4] 

4. GLOBAL CLUSTER TEMPERATURE PROFILE 

The azimuthally averaged gas temperature profile was de- 
rived from the ACIS-S and ACIS-I data by extracting spectra 
in the annular regions indicated in Fig. [3] and is shown in the 
left panel of Fig. [4] The annular bins and temperature mea- 
surements are detailed in Table [2] The outer two temperature 
points represent the farthest temperature measurements from 
the cluster center obtained with Chandra data at present. We 
note that the temperature profile peaks around 180", just in- 
side the inner edge of the shock front. However, as discussed 
in Sect. [6] it is unlikely that this temperature feature is pro- 
duced by the shock as it is consistent with the general shape 
of temperature profiles observed for relaxed galaxy clusters. 

On the other hand, a "plateau " is notable in the tempera- 
ture profile indicating the presence of cool gas in the range 
of radius ^70-150". This peculiar feature has not been noted 
previously in the literature. We thus investigated it in more 
detail thanks to the deep, high-quality Chandra exposure. 

4.1. Evidence for Cool Filaments 

In order to investigate the origin of the plateau seen in the 
temperature profile between ~ 70-150", we made several at- 
tempts to identify regions of cool gas. Our approach was to 
measure a new, "undisturbed" temperature profile by masking 
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TABLE 2 

Temperature Profile derived by Spectral fitting 



Fig. 4. — Left: Azimuthally averaged gas temperature profile derived from the ACIS-S and ACIS-I data by fi tting spectra in the annuli shown 
in Fig. [3] The dotted lines indicate the inner and outer radii of the shock front as determined in Sect. 13.11 Right: Same as left panel, with 
overlaid in red the profile measured after excluding the cool filaments (i.e., after excluding the blue sectors shown in Fig. [5] left panel). 

temperatures by ^0.1-0.2 keV in all bins (except for the very 
central bin, where the temperature is lower by ^0.15 keV) 
up to ^200" . The temperature profile measured by mask- 
ing the cavities is therefore shifted up, maintaining its global 
shape. In particular, the plateau in the temperature profile is 
still present, indicating that the cavities are not a tracer of the 
entire amount of cool gas. We then masked the ~1' bright fil- 
ament stretching from the inner cavity to the center of a larger 
outer cavity in the northeast, which was noted by Nulsen et 
al. (2005). As in the previous attempt, the plateau in the tem- 
perature profile remained, indicating that the known filament 
is not the only repository of the cool gas in Hydra A. 

As a diagnostic of the presence of cool gas we finally used 
the hardness ratio map shown in Fig. [5] which was obtained 
by dividing a smoothed 1.5-7.5 keV image by a smoothed 
0.3-1.5 keV image. The dark regions are indicative of low 
gas temperatures. They agree with the arm-like structures of 
cooler gas extending towards the north and the south in the 
temperature map of Simionescu et al. (2009a). Based on 
a visual inspection of the hardness ratio map, we selected a 
simple combination of sectors reproducing the shape of the 
supposedly cool filaments (see left panel of Fig. 0. We then 
excluded these regions and extracted the spectra in the same 
radial bins as above (see Fig. |5J, and generated the tempera- 
ture profile shown in red in the right panel of Fig. |4](see Table 
12). As evident from the comparison of the two profiles, the 
plateau has largely been removed and the temperature profile 
is typical of cool core clusters (see also Sect. |6). This clearly 
indicates that the masked regions contain cool gas. A detailed 
spectral analysis of the cool filaments is presented in Sect. [5] 

5. SPECTRAL PROPERTIES OF THE COOL GAS 



Bin 


Radius 


kT 


X 2 /dof 




X 2 /dof 


no. 


(") 


(keV) 




(keV) 




1 


0-13 


? qa+0.04 
*" y -0.04 


763/560 


J - 1J -0.07 


453/430 


2 


13-21 


3 07+ 04 


638/590 


3 26+* 06 
J ' zo -0.06 


520/466 


3 


21-29 


3 9Q+0.05 
J -"-0.04 


607/598 


3 34+0.05 
J -"-0.05 


589/566 


4 


29-37 


3 3q+0.05 

^■•"-o.os 


663/600 


3 cc+0.07 
J ' JJ -0.07 


628/565 


5 


37-46 


3 61+ 06 


681/622 


J -' D -0.07 


605/584 


6 


46-58 


3 60 +0 06 
J - -0.06 


644/607 


3 7^+0.07 
J -' J -0.07 


615/578 


7 


58-72 


3 62+ 007 


737/623 


3 o 1+0.07 
J - M -0.07 


618/593 


8 


72-89 


3 c ,+0.07 
J - J -0.07 


656/627 


3 73+0.10 

J -' J -o.io 


512/502 


9 


89-107 


3 c ,+0.07 
JJ -0.07 


748/629 


-> os+o.io 


538/513 


10 


107-126 


3 co+0.07 
J,J -0.07 


698/634 


3 09+O.IO 


507/508 


11 


126-145 


T CO+0.07 

- , ' JO -0.07 


713/648 


-3 7O+0.11 


491/511 


12 


145-168 


3 7 Q+0.07 
J - "-0.07 


744/674 


3 q 1+0 08 

• , - yi -o.m 


672/643 


13 


168-201 


4 Q9+0.08 
-0.08 


719/709 


4 (19+0.08 
^• uz -0.08 


719/709 


14 


201-257 


3 Q 1+0.08 
J - 7 -0.08 


1125/1124 


3 q 1+0 08 
J - ?1 -0.08 


1125/1124 


15 


257-339 


3 70+°' 11 


1252/1177 


3 7n+0.n 


1252/1177 


16 


339-497 


-i co+0.19 
J.38_ 19 


1140/1150 


3 S8+°- 1EI 
J - JO -0.19 


1140/1150 


17 


497-812 


2 96+ 30 
-0.29 


1717/1676 


2 96+0- 30 


1717/1676 



NOTE. — Results of the spectral fitting in concentric 360° annular re- 
gions (shown in Fig. [3j in the 0.5-8.0 keV energy range using the XSPEC 
wabsxapec model. The absorbing column density is fixed to the Galactic 
value (N}i = 4.68 X 10 20 cm~ 2 ), while the temperature (in keV) and metallicity 
(non reported here) are left as free parameters. Error bars are at the 90% confi- 
dence levels on a single parameter of interest. The values of kT ma ^ are measured 
in the same concentric annular regions after excluding the cool filaments (i.e., 
by fitting spectra extracted after masking the blue sectors shown in Fig. \5\ left 
panel). See Sect.[4|for details. 



such regions and then compare it with the global temperature 
profile measured above. We first masked the cavity regions (as 
indicated in Wise et al. 2007), finding systematically higher 



We investigate here the spectral properties of the gas which 
produces the plateau seen in the global cluster temperature 
profile. In particular, we focus on the gas located in the 
range of projected distances from the center of 72-145" (~ 
76-152 kpc), which corresponds to the radial range where the 
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Fig. 5. — Left: High-contrast hardness ratio map obtained by dividing a 1.5-7.5 keV image by a 0.3-1.5 keV image. Each raw image was 
smoothed with a 10 pixel (5") gaussian before the division. The images were extracted from the merged cleaned event files for ObsID's 4969 
and 4970. Regions in black are indicative of low temperature gas. The blue sectors indicate the selected regions of supposedly cool filaments. 
These regions have been excluded in measuring the temperature profile shown in the right panel of Fig. [4] (red points). The red contours 
outlining the 330 MHz radio emission from Lane et al. (2004) are also shown for comparison. Right: Similar to left panel. Overlaid in white 
are the sectors used to study the spectral properties of the cool filaments (located between radii 72-145" ~ 76-152 kpc, see Sect.|5]for details). 



TABLE 3 

Spectral Analysis of the Cool Filaments 
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Sector kT(keV) EM(xl(T 3 ) Z (solar) 


X 2 /dof [xll 


kTi (keV) EMi(xl(T 3 ) kT 2 (keV) EM 2 (xl(T 3 ) Z (solar) 


X 2 /dof F stat 



WNW 


■2 O7+0.10 


1 


60 +0.03 
ou -0.03 


37+0 05 
u ' J '-0.05 


505 


1/521 


[0 


97] 


2 


1 7+0.48 
1 -0.47 


0.73^ 


6 40+ L89 
°-^ u -1.74 


0.94*$ 





90+O.O6 
-0.05 


484 


6/519 


11.0 


NNW 


3 94+O.O6 
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36 +0 04 
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467 
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32.0 


NNE 


3 , 4+0.08 
J ' 1H -0.08 


1 


93+0.03 
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41+0.05 


483 


0/460 


[1 


05] 
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i 9+0.23 
lz, -0.21 
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0M$* 





30+0.05 
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ENE 


3 39+0.08 
■'■•"'-0.08 
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i c+0.02 
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49+O.O6 
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[1 


00] 


2 


94+0.37 
^M>.47 
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0.48^° 





34+0.06 
-^-0.05 
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11.4 
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3 45+0.10 
J, ^ J -0.10 
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3 fi +0.03 


37+0 05 
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, fi +0.86 
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4 49+14.9 


n 09+O.33 
u - yz -0.66 





39+0.05 
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2/481 


5.4 
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3 3 4+0.07 
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NOTE. — Results of IT model and 2T model spectral fitting in the 0.5-8.0 keV energy range. The annular sectors of spectra extraction are indicated in the right panel of Fig. [5] The 
normalizations (EMs) are in XSPEC units of 10~' /7 c "pV/47t[Da(1 + z)] 2 . Errors are at the 90% confidence levels on a single parameter of interest. The last column shows the result of 
theF-test. See Sect.[5]for details. 



plateau in the temperature profile is most evident (i.e., bins 
no. 8-11, see Fig. [4] and Table [2J. We divided the annulus 
from 72-145" into 8 sectors, each having an angular width of 
45°, obtaining the regions labeled as WNW (west-northwest), 
NNW (north-northwest), NNE (north-northeast), ENE (east- 
northeast), ESE (east-southeast), SSE (south-southeast), SSW 
(south-southwest), WSW (west-southwest) in the right panel 
of Fig. [5] We extracted the spectra in these sectors and com- 
pared two different spectral models. The "IT model" is the 
absorbed apec model already used above to derive the global 
temperature profile. The free parameters are the temperature, 



kT, the metallicity, Z, and the normalization (emission mea- 
sure, EM). The "2T model" includes a second thermal emis- 
sion component (apec+apec) and has 2 additional free pa- 
rameters: the temperature, and the normalization, EM2, 
of the second component (the metallicities of the two compo- 
nents are linked). 

The best-fitting parameter values and 90% confidence 
ranges derived from the fits to the annular spectra in sectors 
are summarized in Table [3] Although the improvement of 
adding a second thermal component is formally significant 
according to the F-test, our results show that the quality of 
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Fig. 6. — IT model (left) and 2T model (right) fit to the spectra extracted in the SSE sector. See Table|3]for best-fitting parameter values. 



1.2 









/ 3(7 c.l. " 




r 90% c.l. ■ 







AkT (keV) 

Fig. 7. — Reduced chi square of the IT model fit to mock 2- 
temperature spectra as a function of the separation AT of the 2 ther- 
mal components. The best-fit value of the single temperature fit is 
kT = 3.24 ± 0.06 keV. 

the data is not generally sufficient to demand a model more 
complex than the IT model. In fact, in most sectors, the IT 
model already produces a very good fit (reduced chi squared 
xt ~ 1) and therefore a more complicated model appears un- 
necessary. We also note that the second thermal component is 
poorly constrained, with temperature errors >25% and up to 
300%. Only in the SSE sector is the reduced chi squared of the 
IT model unacceptable at 90% significance, and the statistical 
improvement obtained by introducing an additional emission 
component compared to the single-temperature model is the 
most significant according to the F-test. The improvement 
of the 2T model over the IT model in this sector is also ev- 
ident from the residuals of the fits in Fig. [6] We can there- 



TABLE4 

Spectral Analysis of the Cool Filaments (cont.) 



2T model with kT 2 = 4 keV 
Sector Wi (keV) EM^xlCr 3 ) EM 2 (xl(r 3 ) x 2 /dof F stat 
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1 -'^-0.04 
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NOTE. — Results of 2T model spectral fitting in the 0.5-8.0 keV energy range. The 
temperature value of the second thermal component is fixed to 4 keV. The annular sectors 
of spectra extraction are indicated in the right panel of Fig.[5l The normalizations (EMs) 
are in XSPEC units of I0~ u n c ti p V /4tt[D a (1 + z)] 2 . Errors are at the 90% confidence 
levels on a single parameter of interest. The last column shows the F statistic value 
testing the improvement of the 2T model over the IT model in TablefJ] See Sect.|5]foi- 
details. 

fore conclude that, confirming the hardness ratio map (Fig. |5J 
right panel), we find spectral evidence for multiphase gas in 
the SSE sector with a hot component at 4.92^^ keV and 
a cool component at 1.59^23 keV. Assuming that the two 
spectral phases are in pressure equilibrium in the same vol- 
ume, the ratio of the volumes they occupy is estimated as 
Vi/V 2 = (EM1/EM2) • (kTi/kT 2 ) 2 , so the filling factor of the 
cool gas is ^0.04. 

By contrast, we do not find clear spectral signatures of cool 
gas in the sectors NNW, NNE, and ENE, as expected from 
a visual inspection of the hardness ratio map. However, the 
lack of spectral evidence for multiphase gas could be due 
to the limitations of our data. Indeed, due to the relatively 
limited spectral resolution of Chandra, the detection of two 
different thermal components demands a significant tempera- 
ture separation. The temperature difference required to have 
a marked effect on a single-phase thermal fit at the 90% con- 
fidence level is determined as follows. Using the response 
matrices, background, and numerical information of a real 
spectrum extracted in an arbitrary sector, we simulated spec- 



8 



Gitti et al. 



1.5 



V 



0.5 




1.5 



V 



E- 



0.5 




Fig . 8 . — Left: Temperature profile measured for Hydra A (red triangles, corresponding to the profile shown in the left panel of Fig. |4j overlaid 
on the temperature profiles of a sample of 12 relaxed clusters presented by Vikhlinin et al. (2005). The temperatures are scaled to the cluster 
emission-weighted temperature excluding the central 70 kpc regions. By extracting the global spectrum of Hydra A, after masking the central 
67", we measured a value < Tx >= 3.76 ± 0.03 keV. The profiles for all clusters are projected and scaled in radial units of the virial radius r,„>, 
estimated from the relation r vlI = 2.74Mpc^/< Tx > /10keV (Evrard et al. 1996). Right: Similar to the left panel, but with the temperature 
profile measured for Hydra A after masking the cool filaments (corresponding to the profile shown in the right panel of Fig. |4j. 



tra with two thermal components separated by AT around the 
best-fit value of the single temperature fit to the real spectrum. 
We then fitted a single temperature apec model to the mock 
2-temperature spectrum. We repeated this exercise, increas- 
ing the separation between the two temperature components 
until xt exceeds the 90% confidence range. For comparison, 
we performed this procedure by starting from the real spectra 
in different sectors indicated in Fig. [5] and found consistent 
results. An example of such an analysis is plotted in the left 
panel of Fig. [7] where the dotted and dashed lines are the 90% 
and 3cr limits. We found that a ^2.8 keV (3.0 keV) separation 
is necessary to exclude the presence of single-phase gas at the 
90% (3cr) confidence for our data. The fact that the two spec- 
tral components detected in sector SSE are separated by ^3.3 
keV is consistent with this result. We note that kT\ found by 
the 2T model fit in the SSE sector is the lowest among all of 
the sectors, falling in the energy range which comprises most 
of the counts and thus being more easily detectable. 

Despite this, since a temperature higher than 4 keV is not 
observed at any radius in the cluster (Fig. |4), the second ther- 
mal component found by the 2T model fit appears unrealis- 
tically hot so we performed a new spectral fit with the 2T 
model (apec+apec) keeping the temperature of the second 
thermal component fixed at 4 keV. Such a model has only one 
additional free parameter than the IT model: the normaliza- 
tion, EM2, of the second component (the metallicities of the 
two components are linked). The best-fitting parameter val- 
ues and 90% confidence levels derived from the fits to the 
annular spectra in sectors are summarized in Table |4] The F 
statistics for the improvement over the IT model, shown in 
the last column, indicate where the addition of a second ther- 
mal component is most significant (sectors SSE, NNW, NNE, 
and ENE). In agreement with the hardness ratio map (Fig. |5), 



our spectral analysis therefore supports the presence of mul- 
tiphase gas along the filaments. Interestingly, such cool fil- 
aments follow the morphology of the powerful central radio 
source nicely, although the western part of the southern radio 
lobe appears devoid of cool gas. The SSW sector is indeed the 
hottest and is also the region where the presence of a second 
thermal component is least significant. This sector lies at the 
location where the southern radio lobe appears to fold back on 
itself (Lane et al. 2004). The properties of the cavity are also 
consistent with a sharp bend in the southern jet there (Wise et 
al. 2007). 

Finally, we also attempted to map the emission measure 
distribution by fitting more complicated multi-phase spectral 
models (such as apec+apec+apec+apec, with fixed tem- 
peratures and free normalizations). However, we found our 
data are inadequate due to limited statistics and the spectral 
resolution of Chandra. These detailed spectral studies will, 
hopefully, be possible in the future with the spectral capabili- 
ties of the International X-ray Observatory (IXO). 

6. DISCUSSION 
6.1. Scaled Temperature Profile and Shock Front 

We show in the left panel of Fig. [8] the observed 
temperature profile of Hydra A (red triangles), scaled 
by the virial radius (estimated from the relation r v j, = 
2.74MPCA/K T x > /10keV, Evrard et al. 1996) , overlaid 
on the scaled temperature profiles of a sample of 12 relaxed 
clusters observed with Chandra (Vikhlinin et al. 2005). The 
temperatures are scaled to the emission-weighted cluster tem- 
perature <T X >, measured after excluding the central 70 kpc 
region which is usually affected by radiative cooling. We note 
that the plateau of cool gas, already discussed in Sect. [4] 
stands out from to the typical temperature profile of clusters. 
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TABLE 5 
Deprojection analysis 



Shell kT cool EM cool kT hot EM hot V, ot f coo { « e ,cool M cool / hot n c . hot M hot 

(") (keV) (xl(T 3 ) (keV) (xlO~ 3 ) (cm 3 ) (cirT 3 ) (M Q ) (cm" 3 ) (M Q ) 

Inner l-T^oi 3 - 07 4 0I 1)04 13 3 5.35 x 10™ 0.043 0.032 6.86 x 10 10 0.957 0.014 6.68 xlO 11 

(0-72) 

Middle 1 °g 1.64 4.01^ 10.9 3.85 X 10 71 0.023 0.012 9.90 X 10 10 0.977 0.0047 1.65 XlO 12 
(72-145) 

Outer - 4.01+jj ^ 12.5 5.20 X 10 72 1 0.0136 6.59 X 10 12 

(145-339) 

NOTE. — Results of the deprojection analysis in three concentric 360° annular regions (indicated in the first column) in the 0.5-8.0 keV energy range 
using the XSPEC pro jctxwabsx (apec+apec) model. The normalization of the cool component in the outer shell is fixed to zero, and the tempera- 
ture values of the hot components are linked between all shells. The fit gives x 2 /dof = 4993.5/2811. The normalizations, EM coo i(h i), are in XSPEC units of 
10 _14 n e €O oi(hot)WpV c ool(hoi)/4'?r[DA(l +^)] 2 . Errors are at the 90% confidence levels on a single parameter of interest. The ratio of the relative volumes occupied is 
estimated as Vcool /Vl« = (EM L . ol/EM h oi) ■ (kT cml / kT hol )- , and the filling factor /is defined as /„„i(hoi) = Kool(hoi)/V[oi, with Vj , = V m i + V h oi ■ The masses are thus 
estimated from the relation: M = p/Vjoi, where we use the conversion from electron number density to gas density p = 1.83/i« e m//. 



In particular, lying just inside the typical temperature peak, 
such a plateau emphasizes the temperature peak in Hydra A, 
making it look higher than it actually is. Since this region is 
also where the shock front is located, it makes it difficult to 
distinguish between a temperature jump due to the shock and 
the typical temperature peak of the cluster. 

In fact, once we overlay the scaled temperature profile of 
Hydra A measured after masking the cool filaments (right 
panel of Fig. it is remarkable how well it agrees with the 
general shape of the temperature profiles observed for relaxed 
clusters (see right panel of Fig. |8}. Therefore, it is unlikely 
that the temperature peak observed in Hydra A is produced by 
the shock. In particular, the spectroscopic detection of a tem- 
perature rise in the regions immediately inside of the shock 
front is complicated by the underlying rising temperature pro- 
file of the global cluster atmosphere. In Fig. |9]we show some 
examples of such temperature measurements in various sec- 
tors (see Sect. |3.21 i: the post-shock regions are found to be 
hotter than the corresponding pre-shock ones, although due 
to the large error bars, the pre-shock and post-shock tempera- 
tures are consistent. 

6.2. Evidence for Gas Dredge-Up along the Cool Filaments 

The fact that the scaled temperature profile of Hydra A 
measured after masking the filaments agrees with the general 
shape of the temperature profiles observed for relaxed clus- 
ters (right panel of Fig. [8} is a clear indication that these fil- 
aments are responsible for the plateau of cool gas. In Section 
[5] we performed a detailed (projected) spectral analysis and 
found evidence for multiphase gas, which may have been up- 
lifted with the radio lobes, giving rise to the cooler filaments 
and plateau. In order to estimate the mass of the cool gas 
in the filaments, we must correct for the effects of projection 
along the line-of-sight. Starting from the annuli used to de- 
rive the global temperature profile, we binned them together 
in order to obtain three shells suitable for a simple deprojec- 
tion analysis of the cool filaments. In particular, using three 
annuli: inner (0-72" -0-76 kpc ), middle (72-145" -76-152 
kpc) and outer (145-339" —152-360 kpc), we performed a de- 
projection analysis of the absorbed 2T model with the XSPEC 
pro jet model. Under the assumption of spherical shells of 
emission, this model calculates the geometric weighting fac- 
tors, according to which the emission is redistributed amongst 
the projected annuli. The outer shell was assumed to include 
only gas at the ambient temperature, which is linked to the 
temperature values of the hot components in the inner and 
middle shells. The detailed results of the deprojection analy- 



sis are shown in Table [5] 

We found that the hot phase has a temperature of 4 keV, in 
agreement with the observed global temperature profile (Fig. 
|4), and that the cool component is, on average, at —1.6 keV. 
By taking into account the relative filling factor, /, calculated 
under the assumption that the two thermal components are in 
pressure equilibrium in each shell, in the middle shell we esti- 
mated a mass of cool gas of 9.9 x 10 I0 M Q (and a mass of hot 
gas of 1.7 x 10 12 M Q ). We note that this estimate is consider- 
ably larger than the mass of cool gas reported by Simionescu 
et al. (2009a). The discrepancy between the two measure- 
ments can be ascribed to the fact that we included the mass 
content of the mushroom-cap (namely, sectors ENE, NNE, 
NNW in Fig. [5]). Using the entropy of the cooler phase and 
comparing to the entropy profile of the Hydra A cluster mea- 
sured by David et al. (2001), we could estimate where it orig- 
inated. From the density and temperature values measured in 
the middle shell (see Table |5), the cool gas has an entropy 

S = Wcooiw e cooi ~30 keV cm 2 . As shown by the profile of 
David et al. (2001), gas at such entropy is located around 30 
kpc from the cluster center. Therefore, it is very likely that 
the cool gas, which is now observed in the filaments at —75- 
150 kpc, was lifted from the central 30 kpc. In this case the 
most obvious way to move the gas is by some form of en- 
trainment or dredge up of cool material from the center of the 
cluster in the rising lobes. Churazov et al. (2000) have indeed 
shown that hot buoyant bubbles produced by the expanding 
central radio source can dredge up cool material from the cen- 
ter of the cluster. The mass of cool gas at —75-150 kpc is 9.9 
x 10 10 M Q , which is —60% of the total mass of gas remaining 
within 30 kpc (1.7 xl0 n M Q , David et al. 2001). Assuming 
such a mass was lifted out of the central 30 kpc by a con- 
tinuous outflow or a series of bursts from the nucleus of the 
central galaxy over the past 200-500 Myr (which created the 
cavity system, Wise et al. 2007), it would amount to outflows 
of a few hundred M Q yr" 1 . There would thus be a develop- 
ment of a convectively unstable region that can significantly 
reduce the net inflow of cooling gas (see also David et al. 
2001, Nulsen et al. 2002). Therefore our results show that 
the AGN feedback in Hydra A is acting not only by directly 
re-heating the gas, but also by removing a substantial amount 
of potential fuel for the supermassive black hole (SMBH). 

The energy required to lift the cool gas gives a lower limit 
on the amount of AGN outburst energy deposited in the ICM. 
This value can be estimated by calculating the variation in 
gravitational potential energy during the lifting process. If we 
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FlG. 9. — Temperatures measured in the pre-shock (blue) and post-shock (red) regions of some of the sectors shown in Fig. \\] Radial error bars show the limits of the regions. The 
underlying global cluster profile, measured after masking the cool filaments corresponding to the profile shown in the right panel of Fig. [4] is shown in black. The yellow dashed line 
indicates the radius of the shock front in each sector, while the dotted lines indicate the inner and outer radii of the shock front as determined in Sect. |3.1| The bottom panels are for 
binned sectors. 
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assume that the undisturbed ICM is approximately isothermal 
with sound speed c s w 1000 km s" 1 and is in a hydrostatic 
configuration with density profile p(r), we can calculate this 
quantity as (Reynolds et al. 2008) 



A£ = ^2iiln 

7 



(1) 



where M coo i is the lifted mass, p, and pf are the initial and fi- 
nal densities of the surrounding ICM, and 7=5/3 is the ratio of 
specific heat capacities. From the density profile presented by 
David et al. (2001) we estimated that the energy required to 
lift the cool gas is >2.2 x 10 60 erg. This value is comparable 
to the work required to inflate all of the cavities against the 
surrounding pressure (Wise et al. 2007) and is —25% of the 
total energy of the large-scale shock (Nulsen et al. 2005). Al- 
though we find evidence for a more extended gas dredge-up 
than previously estimated, there is a remarkable correlation 
between the cool filaments studied here and the metal-rich fil- 
aments in the iron-abundance maps measured by Simionescu 
et al. (2009a) and Kirkpatrick et al. (2009). This is consistent 
with the idea that Hydra A's powerful radio source is able to 
lift cool, metal-rich gas from the central region and distribute 
it throughout the X-ray atmosphere of the cluster. A similar 
effect is observed in M87 (Simionescu et al 2008, Werner et 
al. 2010). We finally note that by summing our estimates 
of the integrated mass of cool gas, M coo i(< 152kpc) — 1.7 
xl0 n M Q , and hot gas, M hot (< 152kpc) -2.3 xl0 12 M Q , in 
the inner and middle shells we measure a total mass of gas 
M gas (< 152kpc) —2.5 xl0 12 M Q which is in agreement with 
the gas mass profile measured by David et al. (2001), so our 
general picture is self-consistent. 

7. SUMMARY 

The main results of this work can be summarized as fol- 
lows: 

• We studied the azimuthal properties of the weak (Mach 
number —1.2-1.3), large-scale shock and attempted 
to measure the temperature jump associated with the 



shock in different directions. The errors in the temper- 
ature measurements are too large to constrain the tem- 
perature jump caused by the shock. Furthermore, we 
note that the detection of a temperature rise in the re- 
gions immediately inside of the front is complicated by 
the underlying rising temperature profile of the cluster 
atmosphere. 

• Our detailed spectral analysis indicates the presence of 
multiphase gas along soft filaments seen in the hard- 
ness ratio map. The cooler gas has a significant impact 
on the radial temperature profile of the cluster. After 
correcting for the effect of the cool filaments, Hydra 
A's temperature profile is consistent with the form ob- 
served in relaxed galaxy clusters. Thus it is unlikely 
that the observed temperature peak is produced by the 
shock. 

• The cool filaments are direct evidence for substantial 
dredge-up of low entropy material by the rising lobes: 
—60% of the gas mass remaining in the central 30 kpc 
has been lifted along the cool filaments to the observed 
current position of 75-150 kpc. The outflow amounts 
to a few hundred M Q yr" 1 , which is comparable to the 
inflow rate formerly estimated for the cooling flow. The 
energy required to lift the cool gas is —25% of the total 
energy of the large-scale shock and is comparable to the 
work required to inflate the cavities. 
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